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Progress in extracting excited-state baryon masses in lattice QCD using large sets of spatially-extended operators 
is presented. The use of stochastic estimates of all-to-all quark propagators with variance reduction techniques is 
described. Such techniques are crucial for incorporating multi-hadron operators into the correlation matrices of the 
hadron operators. 



1. INTRODUCTION 

O ; 

Experiments show that many excited-state hadrons exist, and there are significant experimental efforts to map out 
q i the QCD resonance spectrum, such as Hall B and the proposed Hall D at Jefferson Lab, ELSA associated with the 
University of Bonn, COMPASS at CERN, PANDA at GSI, and BESIII in Beijing. Hence, there is a great need for 
,_|_ ■ ab initio determinations of such states in lattice QCD. 

To extract excited-state energies in Monte Carlo calculations, correlation matrices are needed and operators with 
very good overlaps onto the states of interest are crucial. To study a particular state of interest, all states lying 
below that state must first be extracted, and as the pion gets lighter in lattice QCD simulations, more and more 
multi-hadron states will lie below the excited resonances. To reliably extract these multi-hadron states, multi-hadron 
operators made from constituent hadron operators with well-defined relative momenta will most likely be needed, 
and the computation of temporal correlation functions involving such operators will require the use of all-to-all quark 
propagators. The evaluation of disconnected diagrams will ultimately be required. Perhaps most worrisome, most 
excited hadrons are unstable (resonances), so the results obtained for finite-box stationary-state energies in lattice 
k*" ' QCD must be interpreted carefully. 

In this talk, progress by the Hadron Spectrum Collaboration in extracting excited-state baryon masses in lattice 
I/"} QCD using large sets of spatially-extended operators is presented. The use of stochastic estimates of all-to-all 
quark propagators with variance reduction techniques is described. Such techniques are crucial for incorporating 
' multi-hadron operators into correlation matrices. 

O 2. EXCITED STATIONARY STATES IN LATTICE QCD: KEY ISSUES 

Capturing the masses of excited states requires the computation of correlation matrices Cij(t) = (0|T<I>i(t)<I>t(0)|0} 
associated with a large set of N different operators &i(t). It has been shown in Ref. 1] that the N principal effective 
' masses W a (t), defined by 



in 



W a {t) = ln 



^a{t, to) 
X a (t + l,t ) 



where \ a (t, to) are the eigenvalues of C(to)" 1 ^ 2 C(t) C(to)" 1 / 2 and to < t/2 is usually chosen, tend to the eigenenergies 
of the lowest N states with which the N operators overlap as t becomes large. When combined with appropriate 
analysis methods, such variational techniques are a particularly powerful tool for investigating excitation spectra. 

The use of operators whose correlation functions C{t) attain their asymptotic form as quickly as possible is crucial 
for reliably extracting excited hadron masses. An important ingredient in constructing such hadron operators is 
the use of smeared fields. Operators constructed from smeared fields have dramatically reduced mixings with the 
high frequency modes of the theory. Both link-smearing and quark-field smearing should be applied. Since excited 
hadrons are expected to be large objects, the use of spatially extended operators is another key ingredient in the 
operator design and implementation. Fig. [T] shows the different spatial configurations we use, which effectively build 
up the necessary orbital and radial structures of the hadron excitations. The basic building blocks in all of our hadron 
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Figure 1: The spatial arrangements of the extended three-quark baryon and quark-antiquark meson operators. Solid circles 
are smeared quark fields, open circles in the mesons are smeared "barred" antiquark fields, open circles in the baryons indicate 
a Levi-Civita color coupling, and solid line segments indicate covariant displacements. 



operators are covariantly-displaced quark or antiquark fields. These are first combined to have the appropriate flavor 
structure and color structure, then group-theoretical projections are applied to obtain operators which transform 
irreducibly under all lattice rotation and reflection symmetries. A more detailed discussion of these issues can be 
found in Ref. 0. 

A first glimpse of the higher-lying nucleon spectrum in lattice QCD was provided by the Hadron Spectrum Col- 
laboration in Ref. These first results, shown in Fig. [2j were on small 12 3 x 48 anisotropic quenched lattices with 
a very heavy pion. Results for both the nucleons and A-resonances on 239 quenched configurations on a 16 3 x 64 
lattice and 167 quenched configurations on a 24 3 x 64 lattice using an anisotropic Wilson action with spatial spacing 
a s ~ 0.1 fm, a s /a t ~ 3, and a pion mass ~ 490 MeV appeared during the past yearjj). These masses have also 
been determined in the past year using 430 Nf — 2 configurations on a 24 3 x 64 lattice with a stout-smeared clover 
fermion action and a Symanzik-improved anisotropic gauge action Q. The results for a pion mass = 400 MeV, 
spacing a s ~ 0.1 fm and a s /a t ~ 3 are shown in Fig. [2] The low- lying odd-parity band shows the exact number of 
states in each channel as expected from experiment. The two figures show the splittings in the band increasing as 
the quark mass is decreased. At these heavy pion masses, the first excited state in the G\ g channel is significantly 
higher than the experimentally measured Roper resonance. It remains to be seen whether or not this level will drop 
down as the pion mass is further decreased. Most of the levels in the right-hand plot lie very close to two-particle 
thresholds. The use of two-hadron operators will be needed to go to lighter pion masses. 



3. MANY-TO-MANY QUARK PROPAGATORS 

To study a particular eigenstate of interest, all eigenstates lying below that state must first be extracted, and 
as the pion gets lighter in lattice QCD simulations, more and more multi-hadron states will lie below the excited 
resonances. The correlation functions of such operators require estimates of the quark propagators from all spatial 
sites on a time slice to all spatial sites on another time slice. Computing all such elements of the propagators exactly 
is not possible (except on very small lattices). Some way of stochastically estimating them is needed. 

Random noise vectors r\ whose expectations satisfy E(r]i) = and E(r)iT]*) = Sij are useful for stochastically 
estimating the inverse of a large matrix M as follows. Assume that for each of Nr noise vectors, we can solve the 
following linear system of equations: Ml' r ' = rf r ' for X^ r \ Then X^ = M^ 1 rf r \ and 

E{X iV *) = M7 k \ kV *) = £ M^EfrkV*) = M^Sv = My 1 . (1) 

k k k 

The expectation value on the left-hand can be approximated using the Monte Carlo method. Hence, a Monte Carlo 
estimate of M~^ 1 is given by 

, N R 

M^ra lim — V AfV W *, where MX M = r/ r >. (2) 

Nr—>oo N a — ' J 
r— 1 

This equation usually produces stochastic estimates with variances which are much too large to be useful. 

Progress is only possible if stochastic estimates of the quark propagators with reduced variances can be made. A 
technique of diluting the noise vectors has been developed which accomplishes such a variance reduction A given 
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Figure 2: (Left) Nucleon spectrum from 200 quenched configurations on a 12 3 x 48 anisotropic lattice using the Wilson gauge 
and quark actions with a a ~ 0.1 fm, a 3 /at ~ 3.0 and m w ~ 700 MeV from Ref. [3J. (Right) Nucleon spectrum from 430 
iV/ = 2 configurations on a 24 3 x 64 lattice using a stout-smeared clover fermion action and Symanzik-improved gauge action 
with a s ~ 0.1 fm, a s /a t ~ 3, and mi = 400 MeV from Ref. 0. 



dilution scheme can be viewed as the application of a complete set of projection operators. To see how dilution works, 
consider a general N x N matrix M having matrix elements My. Define some complete set of N x N projection 
matrices P<» which satisfy p( a )pW = 6 ab P^ with J2 a p{a) = 1 and p(a)t = p(a) ■ Thcn observe that 

a a a 

= E M * P wE(v*4)P$ Y, M ik E ( P tU>r,*>P$) ■ (3) 

a a 

Define 77^ = P^Vk' and ry^* = f]pP^ — Pjf*Vj'i an d further define At a J as the solution of M ik X^ = t]^\ then 
we have 

M J = E ^foN" 1 *) = J2 E( X \ a W; ] *)- (4) 

a a 

Although the expected value of ^ a r^rj^* is the same as ??fc?y*, the variance of ^ a fjj^nj-^* is significantly smaller 
than that of r\k"tfj- For both Z4 and t/(l) noise, we have Var(Re(riiri*)) = Var(Im(r]iT]j)) — |(1 — <5y). Although 
the variance is zero for i = j, there is a significant variance for all i j. The dilution projections ensure exact zeros 
for many of the off-diagonal elements, instead of values that are only statistically zero. In other words, many of the 
i =/= j elements become exactly zero. 

Of course, the effectiveness of the variance reduction depends on the projectors chosen. A particularly important 
dilution scheme for measuring temporal correlations in hadronic quantities is "time dilution" where the noise vector 
is broken up into pieces which only have support on a single time slice. Spin and color dilution are two other easy- 
to-implement schemes, and various spatial dilution schemes are possible. These various dilution projectors can also 
be combined to make hybrid schemes. 

A comparison of the different dilution schemes is shown in Fig. [3J These computations are dominated by the 
inversions of the Dirac matrix, so using the number of matrix inversions A; nv to compare computational efforts is 
reasonably fair. The advantage in using increased dilutions compared to an increased number of noise vectors with 
only time dilution is evident in the plots. However, this advantage quickly diminishes after time + even/odd-space 
dilution, or time+color, or time+spin dilution. These encouraging results demonstrate that the inclusion of good 
multi-hadron operators will certainly be possible using stochastic all-to-all quark propagators with diluted-source 
variance reduction. 
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C(t=5) for triply-displaced-T nucleon 
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Figure 3: (Left) The relative errors in the correlation function of a single-site nucleon operator for temporal separation t = 5a t 
evaluated using stochastically-estimated quark propagators with different dilution schemes against 1/ivV 2 , where N- lnv is the 
number of Dirac matrix inversions required. The open circle shows the point-to-all error, and the horizontal dashed line shows 
the gauge-noise limit. The black (red) dashed-dotted line shows the decrease in error expected by simply increasing the number 
of noise vectors, starting from the time (time + even/odd-space) dilution point. (Right) Same as the left plot, except for a 
triply-displaced-T nucleon operator. These results used 100 quenched configurations on a 12 3 x 48 lattice. 



4. SUMMARY AND OUTLOOK 



This talk discussed the key issues and challenges in exploring excited hadrons in lattice QCD. The importance 
of multi-hadron operators and the need for all-to-all quark propagators were emphasized. The technology needed 
to extract excited stationary-state energies, including operator design and held smearing, was detailed. Efforts in 
variance reduction of stochastically-estimated all-to-all quark propagators using source dilutions was described. 

Given the major experimental efforts to map out the QCD resonance spectrum, such as Hall B and the proposed 
Hall D at Jefferson Lab, ELSA, COMPASS, PANDA, and BESIII, there is a great need for ab initio determinations 
of such states in lattice QCD. The exploration of excited hadrons in lattice QCD is well underway. 
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